Non-equilibrium relaxation of an elastic string in a random potential 
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We study the non-equilibrium motion of an elastic string in a two dimensional pinning landscape 
using Langevin dynamics simulations. The relaxation of a hne, initially flat, is characterized by 
a growing length, L{t), separating the equilibrated short length scales from the flat long distance 
geometry that keep memory of the initial condition. We show that, in the long time limit, L{t) has 
a non-algebraic growth with a universal distribution function. The distribution function of waiting 
times is also calculated, and related to the previous distribution. The barrier distribution is narrow 
enough to justify arguments based on scaling of the typical barrier. 



The physics of disordered clastic systems has been the 
focus of intense activities both on the theoretical and 
experimental side. Indeed it is relevant in a large number 
of experimental situations ranging from periodic systems 
such as vortex lattices ,1] , charge density waves 1^ ^and 
Wigner crystals Jalto domain walls in magnetic |j, |M 
or ferroelectric 013 systems, contact lines and fluid 
invasion in porous media |lfll | . Because of the competition 
between disorder and elasticity, glassy properties arise, 
and one of the most challenging question is to understand 
their consequences on the dynamics of the system 

Since the system must move by thermal activation 
over the barriers separating metastable states, the steady 
state response to a small external force is a way to probe 
its glassy nature. The glassiness leads to divergent barri- 
ers and thus to a slow response known as creep 12, 13J. 
Experiments 0, 0, 13 as well as microscopic calcula- 
tions of the response have confirmed this creep 
behavior, although questions remain in low dimensions 
about the value of the creep exponent 16]. Much less 
is known about the glassy effects in the case of non- 
stationary relaxation towards equilibrium. Understand- 
ing such non-stationary physics is clearly crucial since 
it gives complementary information on the barriers and, 
for experiments, is needed to describe the many systems 
that are quenched in the glassy state (e.g. by changing 
rapidly the temperature), and have then to relax. The- 
oretical attempts to tackle this problem have been made 
mean field and renormalization group approaches 
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using mean 

[itI 13 Q • Direct application of these results to one 
dimensional domain walls is however difficult. Numeri- 
cal studies, that would give more direct information in 
low dimension, are also difficult since they have to deal 
with ultra long time scales dynamics. Simulations have 
thus been mostly restricted so far to 2-dimensional ran- 
dom Isin g rn odels or 2-dimensional periodic elastic sys- 
tems |2l[ |23| ■ The relaxation of a directed polymer 
has been investigated by local Monte Carlo 

dynamics l26l| , but a precise study of the connection be- 
tween relaxation and the static glassy properties is still 
lacking. 

In this paper we thus study the slow non-equilibrium 



FIG. 1: (a) Typical configurations of the string for different 
times, at T = 0.5. (b) Structure factor for different times at 
T = 0.5, averaged over 1000 disorder realizations. The dashed 
line corresponds to the thermal equilibrium solution which it 
is reached at very long times. 



relaxation of an elastic string moving in a two dimen- 
sional random media. We prepare the string in a flat con- 
figuration and let it relax. We show that the relaxation 
is governed by a characteristic growing length, L{t), sep- 
arating the equilibrated short length scales from the flat 
long distance ones that keep memory of the initial con- 
dition. In the long time limit, L{t) has a non-algebraic 
growth with a universal distribution function. We com- 
pute the distribution of waiting times and thus of barri- 
ers. This later distribution is found to be narrow enough 
to justify the scaling for L{t) based on a typical barrier. 

We consider a string described by a single valued 
function measuring its transverse displacement 
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u from the z axis at time t. The initial condition is flat 
u{z,t = 0) = 0, and we monitor the relaxation towards 
equilibrium. The string obeys the equation of motion: 

jdtU(^z, t) = cd'^u{z, t) + Fp{u, z) + 7](z, t) (1) 

where 7 is the friction coefficient and c the elastic con- 
stant. The pinning force Fp{u,z) = — z) de- 
rives from the random bond disorder potential U (u, z) 
and the thermal noise ri(z,t) satisfies {ri{z,t)) — 
and {r]{z,t)r]{z',t')) = 2-fTS{t - t')S{z - z') where 
(. . .) is the thermal average. The sample to sam- 
ple fluctuations of the random potential are given by 
[U{u,z) - U{u',z')]'^ = -2<5(z - z')R^{u - u') where 
denotes an average over disorder realizations. In the ran- 
dom bond case the correlator R(u) is short ranged. 

To solve numerically we use the method of [l6| . 
We discretize the string along the z direction, z ^ j — 
0, . . . , L — 1, keeping Uj{t) as a continuous variable. A 
second order stochastic Runge-Kutta method is used to 
integrate the resulting equation. To model a continuous 
random potential we generate, for each j, a cubic spline 
U{uj,j) passing through re gula rly spaced uncorrelated 
Gaussian random points |l6l |27| . To characterize the ge- 
ometry of the line during the relaxation we introduce the 
structure factor S{q,t) = s{q,t) = j;{UqUq) where Uq — 

J2j=o '"^(0^"*''"' ^^'^ 1 = 27rn/L with n = 1, . . . , L — 1. 

In the absence of disorder the relaxation of the string 
can be solved analytically and S{q, t) is 

Spureiq.t) ^ S;i^M)[l - exp(-2cg2t/7)] (2) 

where Spl^^^{q) — T/cq"^ is the structure factor at equi- 
librium. From ^ we can separate two regimes: (i) at 
large q the line is equilibrated with the thermal bath and 
its geometry is described by the equilibrium roughness 
exponent ( — 1/2. This behavior can be extracted from 
the power law decay of the structure factor, (ii) 

at small q, however, the string has still a memory of the 
flat initial condition, and the structure factor reaches a 
plateau: Spure{q — > 0+,t) = {2T/j)t. The crossover be- 
tween these two regimes is driven by a unique growing 
characteristic length scale L{t) that can be defined from 
the intersection point of the two limiting behaviors. In 
the pure case Lp^reit) = 2Tr^2ct/j and its power law 
growth defines the dynamical exponent z, as L{t) '-^ t^^^. 

We now discuss our numerical results for the system 
with disorder. We simulate lines of size L = 256, 512, 
1024 with c = 7 = 1. We take i?(0) = 1 and tempera- 
tures ranging from T = 0.1 to T = 0.7. In Fig. ^a) we 
show the typical relaxation of a string. Note that in the 
pure case, for the same parameters, the equilibration of 
a line of size L = 256 occurs after a time t ~ 10^ (see 
Fig. 121). The presence of barriers in the disordered case 
makes the dynamics much more slow, and equilibrium is 
not yet reached at time t = 10^. We show in Fig. ^b) 
the evolution of S{q,t). As in the pure case two regimes 
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FIG. 2: Growing characteristic length scale L{t) of a string of 
size L — 256. (o) symbols correspond to the relaxation of the 
clean system and the dotted line to the analytical result. (•) 
symbols correspond to the disordered case. The continuum 
line is a fit to equation |1] the dashed line is a fit to the power 
law growth at intermediate scales. For 10'' < t < 10® we get 
9 fa 0.49. Inset: exponent 8 extracted from the fit to equation 
|l]in the time interval ti < t < 10®. (A) symbols correspond 
to a system size L = 512, and (o) symbols to L = 256. 

are observed. At short length scales the line has reached 
equilibrium in the random environment and it is charac- 
terized by the well known roughness exponent = 2/3 
|28j . At large length scales a plateau is still present and 
a crossover growing length L{t) can be defined. Quite 
generally the scaling form of S{q, t) can be written as, 

S{q,t) = S'''>{q)G{qL{t)) (3) 

where G{x 0) ~ x^^"^^ and G{x 00) = 1. The ana- 
lytical calculation of L{t) is clearly a non trivial task, but 
a simple estimate can be done relying on phenomenolog- 
ical scaling arguments, based on creep. At low tempera- 
tures the relaxation is dominated by the energy barriers 
U{L) that must be overcomed in order to equilibrate the 
system up to a length scale L. Using the Arrhenius ther- 
mal activation law we can thus express the relaxation 
time t{L) ~ exp[/3?7(L)]. Even if the exact numerical 
determination of U{L) is an NP-complete problem it is 
usually conjectured that the typical barriers of the energy 
landscape scale, asymptotically with L, the same way as 
the free energy fluctuations: U[L) ~ L^, with 6* = 1/3 
for a line. Numerical calculations |2^ and FRG calcu- 
lations 01 seem to confirm this conjecture. Following 
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these arguments we infer that |3C 



L{t) « Lc 



(4) 



where Lc is a characteristic length which can be identi- 
fied with the Larkin length [sj , Uc the associated energy 
scale Uc = U (L^) and to a microscopic time scale. An al- 
ternative form of L{t) would be the power law scaling of 
the clean system, L(t) '~ t^^^, but with a new exponent 
z > 2 taking into account the effect of the energy barri- 
ers. Note that this proposal corresponds to thermally ac- 
tivated motion over barriers scaling logarithmically with 
the size L. Such behavior has been observed in vari- 
ous 2-dimensional disordered systems including periodic 
elastic systems in the so called "marginal glass phase" 
Ellis El Hi- For this model it is possible to show 
that the dynamical exponent takes the form z{T) oc 1/T. 
Moreover the relaxation towards a steady state of an 
elastic string just above the depinning threshold shows 
the same power law behavior with a dynamical exponent 
z < 2 [13. 

We now compare our results with the above different 
scenarios. The growing length scale L{t) can be deter- 
mined from the average structure factor S{q, t) shown in 
Fig. n^b). In practice we define L{t) as the intersection 
between the plateau 5*4 = S{q 0"'",t) and the equili- 
brated structure factor S°'^ q""^^^. The result is shown 
in Fig. 12 Note that the whole time dependence of L{t) is 
described neither by 10} neither by a pure power law. The 
latter scaling can only approximately fit the short time 
relaxation: the fitted dynamical exponent z strongly de- 
creases with increasing temperature and ranges from 20 
to 4. However, for long times, this powerlaw scaling can 
be ruled out due to the observed bending in the log-log 
scale. To be sure that this bending is not an artifact of 
the proximity of the finite size equilibration we verified its 
presence for bigger systems up to a size L — 1024 where 
L{t) <^ L for all considered times. For this reason the 
logarithmic growth seems to be more adequate for long 
times. A two-parameters fit to (0} gives an exponent 9 
which, at long times, becomes size and time independent, 
as shown in the inset of Fig. El 

Although the logarithmic growth law describes well our 
data at long times, we find an exponent 9 « 0.49, bigger 
than the expected value 1/3. If we assume that the dy- 
namics of relaxation is governed by Arrhenius activation, 
this result indicates either a violation of the expected 
scaling of barriers or the presence of non-negligible sub- 
leading corrections in this scaling at the length scales 
spanned by L{t) in our simulations. The inset of Fig. |3 
shows, for different time-windows, the exponent 9. The 
saturation of 9 excludes strong sub-leading corrections 
at least for the largest times reached in the simulations. 
However, the adequacy of the fit with 9 « 0.49 in the 
last three decades is still not enough to exclude the 



presence of logarithmic corrections in the leading term: 
U{L) ^ L^/^log^(L). The latter scenario is consistent 
with the upper bound scaling found numerically in |2^ 
for the barriers separating metastable states of a directed 
polymer in 2-dimensional random media. Such a scaling 
has been shown to also fit well the Monte Carlo re- 
laxation data for a directed polymer. 

The scaling of the barriers U{L) and the subsequent 
evolution of L{t) refer to typical values of U and L{t). 
On the other hand, for broad enough distributions typi- 
cal and mean values can be very different js^ • Therefore, 
the deviations of the numerical data from the predicted 
behavior might be produced by a broad distribution 
of barriers. To check for such a possibility, and to ex- 
tract the barrier distribution, we study the sample to 
sample fluctuations of the various observables. A con- 
venient quantity to compute for each evolving sample is 
the instantaneous value of the structure factor plateau 
s = s{q — > 0+,t) wich is directly related to the growing 
length I ~ ^1/(1+20 As raw data directly confirms, this 
quantity is a stochastic process growing monotonically 
with the time t. Thus, its sample to sample fluctuations 
can be directly related to the distribution of relaxation 
times r and to the statistics of barriers u by assuming 
Arrhenius activation, u log(r) 35]. One obtains |3^ 



(5) 



where u (resp. s) is the sample dependent barrier (resp. 
structure factor plateau) and ^siu) (<I>,((s)) its cumula- 
tive distribution function for a given value of s (resp. u) 
(i.e., <i>s(uo) is the probability to find a barrier u smaller 
than uq, given a fixed value s for the plateau). Fig.|3Ja) 
shows ^u{s) as a function of s for different values of u. 
For all u the distributions are narrow and, on a logarith- 
mic scale, appear just shifted. This suggest the simple 
rescaling s/St, which collapses all the curves as shown 
in the inset. Strikingly, we find that this rescaled func- 
tion for the fully disordered system is indistinguishable, 
at the resolution of our numerical study, from the one 
($„(a; = s / St) = 1 — exp(— x)) of the clean system, and 
from the identical one would obtain for the Larkin model 
[3^ of disorder (despite the fact that this model does 
not have pinning and metastable states). This scaling 
form implies that the sample to sample fluctuations of 
the growing length, l{t), are given by ^u{x = l/Lt) — 
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ith ( 



,i-i-2C 



exp[-aa;^+^'^], with a = r(l-|- jq^) The statis- 
tics of barriers is obtained from Q using the evolution 
St vs u \og{t) of Fig. 12 In Fig. ^h) we show that 
the cumulative distribution <&„(s) derived using the lat- 
ter method indeed coincides with the one obtained from 
a direct analysis of the raw data of s vs u for each sam- 
ple. As for the sample dependent plateau s (for given 
values of u), the distributions of u for given values of 
s are found to be exponentially narrow. Scaling argu- 
ments based on typical values are therefore justified, since 
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FIG. 3: (a) Cumulative distribution function of the plateau 
value s for fixed times t ranging from 5 to 10^. (o) symbols 
indicate the mean value St- Inset: collapse of the cumulative 
distributions in the rescaled variable s/St- (b) Cumulative 
distribution function of the barriers u for s = 10, 20, 40, 80. 
Circles (•) indicate the mean value Ua- Step-lines are obtained 
from raw data, while solid lines are obtained from the rescaled 
cumulative distribution of s, JSJ. 



they can be safely translated directly to the mean values. 
This indicates that the effect of sample to sample fluc- 
tuations cannot explain the deviations of the numerical 
data with respect to the phenomenological predictions 
observed in Fig. [21 and that such deviations must come 
from the scaling of the barriers. Note also that, as visible 
in Fig. |2Ib), the barrier distribution $s(u), contrarily to 
the distribution of plateaux <i>„(s), does not scale with 
u/Us, where Us = (u) is the mean value. Such a scal- 
ing would only work if a pure power law scaling of the 
barriers with length were perfectly verified. The complex 
behavior of <i>s (u) clearly comes from the existence of two 
regimes in the scaling of the barriers as a function of time 
(length) as shown in Fig. [51 However, it remains to be 
understood why the presence of these two regimes does 
not affect the perfect collapse for $u(s) as a function of 
s/St, ranging from the shortest to the longest times. No 
analytical explanation of this fact, nor of the form of the 
corresponding scaling function exists so far. 

At long time, approximate power-law scaling for the 
barriers is recovered, and thus the distribution of barri- 
ers would scale with u/Us- In this case (0 shows that 
the universal function for all the cumulative distribution 
functions would be a stretched exponential. This 



form is different from the one that extremal statistics 
arguments would suggest [s^, prompting for a reexam- 
ination of the physical understanding of the barrier dis- 
tribution in such disordered systems. 
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